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Abstract 

This paper describes a self-contained, automated methodology for active flow control which 
couples the time-dependent Navier-Stokes system with an adjoint Navier-Stokes system 
and optimality conditions from which optimal states, i.e., unsteady flow fields and controls 
(e.g., actuators), may be determined. The problem of boundary layer instability suppres- 
sion through wave cancellation is used as the initial validation case to test the methodology. 
Here, the objective of control is to match the stress vector along a portion of the boundary 
to a given vector; instability suppression is achieved by choosing the given vector to be that 
of a steady base flow. Control is effected through the injection or suction of fluid through 
a single orifice on the boundary. The results demonstrate that instability suppression can 
be achieved without any a priori knowledge of the disturbance, which is significant because 
other control techniques have required some knowledge of the flow unsteadiness such as 
frequencies, instability type, etc. The present methodology has been extended to three 
dimensions and may potentially be applied to separation control, re-laminarization, and 
turbulence control applications using one to many sensors and actuators. 
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1. Introduction 


In the last decade, increased attention has been devoted to the development of tech- 
niques capable of enhancing our ability to control the unsteady flow in a wide variety of 
configurations such as engine inlets and nozzles, combustors, automobiles, aircraft, and 
marine vehicles. Controlling the flow in these configurations can lead to greatly improved 
efficiency and performance, while decreasing the noise levels generally associated with the 
otherwise unattended unsteady flow. Depending on the desired result, one might wish to 
delay or accelerate transition, reduce drag or enhance mixing. There might be a need to 
postpone flow-separation, increase lift or manipulate a turbulence field. Gad-el-Hak 1 and 
Gad-el-Hak and Bushnell 2 provide an excellent introduction to and overview of various 
control methodologies. 

Small improvements in system performance often lead to large payoffs. For exam- 
ple, Butter 3 estimates that a 5 percent improvement in landing maximum lift coefficient 
(Cj(max)) can translate to a 25 percent increase in payload. Cousteix 4 notes that 45 per- 
cent of the drag for a commercial transport transonic aircraft is due to skin friction on the 
wings, fuselage, fin, etc., and that a 10 to 15 percent reduction of the total drag can be 
expected by laminarizing the flow over the wings and the fin. This translates into a reduc- 
tion in fuel requirements, improved performance, and/or increased payload. Muirhead 5 
has shown with a wind tunnel investigation that control of flow separation on a tractor- 
trailor truck can reduce the drag by 30 to 40 percent of the baseline truck configuration. 
This translates into a savings of millions of barrels of fuel per year. 

Encouraged by the potential for huge rewards with what may be a modest input, 
research into ways of achieving the above gains is ever important. In many applications, 
the flow starts from a smooth laminar state, which is inherently unstable, and develops 
instability waves. These instability waves grow exponentially, nonlinearly interacting, and 
lead ultimately to fully developed turbulence or flow seperation. Therefore, one goal of 
a good control system is to inhibit, if not eliminate, these instabilities which lead to the 
deviation from laminar to turbulent flow state. 

1.1 The Wave-Cancellation Concept 

The simplest form of control which might achieve this objective is the wave- 
cancellation approach, based on the premise that the instability mechanisms in low-speed 
transition are dominated by a single instability wave; therefore, cancelling this wave will 
preclude the nonlinear interactions leading to laminar-turbulent transition. The wave- 
cancellation method further assumes that a wavelike disturbance can be cancelled by in- 
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traducing another wave equal in amplitude but opposite in phase, and thus it is mostly 
applicable to systems governed by linear or quasi-linear equations. The key is to deter- 
mine the parameters of the downstream wave which counter (cancel) the evolution of the 
upstream generated wave. Because there are a number of experiments and numerical sim- 
ulations which validate this approach, the wave-cancellation problem is an excellent test 
problem for the new flow-control methodology described in this paper. 

Most of the experiments, aimed at verifying the wave-cancellation concept, were con- 
ducted on the flat plate, except those of Ladd and Hendricks 6 , Pupator and Saric 7 , and 
Ladd, 8 who considered axisymmetric bodies. Many of these experiments were conducted 
in water tunnels. Vibrating wires, 9 hot strips, 10,11 suction and blowing, 7,8 electromagnetic 
generators, 12 and adaptive heating 6 are some of the methods that were used to generate the 
controlling wave. All these input mechanisms give the necessary control of the phase and 
amplitude of the input wave. Among the more successful studies, Milling 9 and Thomas 12 
achieved at least an 80 percent reduction in the input amplitude of the 2-D wave (with 
0.6%-l% amplitude). However, it was not possible to achieve relaminarization, proba- 
bly because of the three-dimensionality of the flow resulting from the interaction between 
background disturbances and the primary 2-D wave. As expected, the studies conducted 
on axisymmetric bodies produced relatively less wave cancellation because these flows are 
inherently three-dimensional. Furthermore, good wave cancellation requires a linear sys- 
tem with constant coefficients. This requirement is clearly violated for flow over a body 
with curvature. It should also be noted that Liepmann and Nosenchuck 11 compared the 
effects of steady and non-steady heating, and found that steady heating demands a 2000% 
increase in energy over an unsteady wave-cancellation technique. Hence, unsteady control 
may be more efficient for flow control applications. 

In addition to the aforementioned experiments, several theoretical (i.e., linear compu- 
tations and theory) and computational studies (i.e, nonlinear simulations) have focused on 
understanding the physics of this wave-cancellation process. Maestrello and Ting 13 pro- 
vided a linear asymptotic analysis to demonstrate the relationship between the input of 
localized disturbances and their effect on the Tollmien-Schlichting (TS) instability waves 
present in the wall-bounded shear flow. They showed that small amounts of local periodic 
heating could excite disturbances that actively control the TS waves which propagate in 
a boundary layer on a flat plate. Analogous to the experiments, several wave input mech- 
anisms were considered. In one of the early Navier-Stokes simulations of active control, 
Biringen 14 used suction and blowing at the wall in a channel flow. He observed approx- 
imately a 50 percent reduction in the amplitudes of the 2-D instabilities and a decrease 
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in the growth of the 3-D instabilities. The Reynolds stress originally generated by the 
waves was all but removed. On the other hand, Metcalfe et al. 15 studied the effect of a 
moving wall on unstable waves traveling in a laminar flow on a flat plate. The simulations 
were based on the Navier-Stokes equations solved within a temporal framework. An energy 
analysis revealed that the wall motion causes the Reynolds-stress term to become negative, 
which implies a feed of energy from the perturbed flow back into the mean flow. In effect, 
this energy analysis showed how a perturbation to an unstable flow can be stabilizing. 
However, an instability wave eventually formed downstream of the control, with the same 
growth rate as the uncontrolled wave. This is a clear indication that the cancellation was 
not complete. Although intuitively obvious, until the work of Bower et al. 16 and Pal et 
ah, 17 it was not known that perfect cancellation could be obtained within the context of 
linear theory (for which the mean flow is independent of the propagating direction). They 
used the 2-D Orr-Sommerfeld equation to study and control instability-wave growth by 
superposition, and showed, within the limits of linear stability theory and the parallel-flow 
assumption, that single and multifrequency waves can be cancelled. The basic conclusions 
drawn by the early experimentalists were confirmed by the studies of Laurien and Kleiser 18 
and Krai and Fasel. 19 They showed that it was possible to delay/accelerate transition by 
superposing disturbances out of/in phase with the primary TS wave. Similar results were 
also reported by Danabasoglu et al. 20 Finally, Joslin et al. 21 performed a numerical ex- 
periment which served to unequivocally demonstrate the link between linear superposition 
and instability suppression. To ensure that linear superposition of individual instabilities 
was, in fact, responsible for the results found in previous experiments and computations, 
they carried out three simulations with i) only the disturbance; ii) only the control; and 
iii) using both disturbance and control, which is the wave-cancellation case. By discretely 
summing the control-only and forcing-only numerical results, they found that this linear 
superposed solution is identical to the wave-cancellation results. These tests deary verify 
the hypothesis that linear superposition is the reason for the previous experimental and 
computational results. 

^From the above experiments, linear computations and nonlinear simulations, several 
common features emerge: (1) It is impossible to achieve perfect wave cancellation unless the 
system is linear, with constant coefficients; (2) the efficiency of wave cancellation decreases 
as the system becomes more nonlinear; (3) as the geometry of the configuration becomes 
more complex, cancellation becomes more difficult; and (4) the current approaches require 
foreknowledge of the instability wave characteristics, such as its frequency and amplitude 
before wave suppression can be realized. 
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1.2 Optimal Control Theory 


The optimal control theory provides an approach which does not require a priori 
knowledge of the flow characteristics. The goal of optimal control theory is to minimize or 
maximize an objective function in a robust manner. When the flow is time-dependent, and 
a strong function of initial conditions, it becomes difficult to establish the precise controls 
that will achieve the desired effect. Wave-cancellation, as discussed above, only works 
well when the input wave has a dominant frequency, and its properties are known. Then 
(either in a passive fashion, or through a feedback mechanism), one seeks to cancel its 
effect while still in a linear regime. In practice, there are many waves, which can interact 
nonlinearly in ways not always known in advance. Rather than try to cancel the incoming 
waves, one seeks appropriate controls in other ways. One means of achieving this, without 
an extensive search over the space of possible controls, is to postulate a family of desired 
controls (e.g., an arbitrary time-dependent amplitude and a specified spatial distribution), 
and an objective function (i.e., stress over a region of the plate). Then, through a formal 
minimization process, one derives a set of differential equations, and their adjoints, whose 
solution produces the optimal actuator profile (among the specified set). While the solution 
to this set of equations cannot be accomplished in real time, the results can be applied 
using standard passive or active control mechanisms. The advantage of this approach is 
that entire collections of controls can be studied simultaneously, rather than one at a time. 

Optimal control methodologies have been recently applied to a variety of problems in- 
volving drag reduction, flow and temperature matching, etc. to provide more sophisticated 
flow control strategies in engineering applications. Computational fluid dynamics (CFD) 
algorithms have reached a sufficiently high level of maturity, generality, and efficiency so 
that it is now feasible to implement sophisticated flow optimization methods, which lead 
to a large number of coupled partial differential equations. Optimal control theory is quite 
mathematical, and its formal nature is amenable to the derivation of mathematical the- 
orems related to existence of solutions and well-posedness of the problem. (Only partial 
results of this type are possible in three-dimensions since, in this case, the Navier-Stokes 
equations themselves do not enjoy a full theoretical foundation; in two-dimensions, a com- 
plete theory is available.) Two recent surveys of the mathematical theories of optimal flow 
control are Gunzburger 22 and Borggaard et al . 23 A mathematical study of a simplified 
problem related to the one considered in this paper can be found in Fursikov et al . 24 

Flows may also be controlled through controllers whose actuation is determined by 
feedback laws. For unsteady flows, including turbulent flows, feedback controls are con- 
sidered, for example, in Choi et al . 25 and Bewley et al. 26 , as well as in Joslin et al . 28 
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Determining optimal feedback laws is a very difficult proposition, especially in the context 
of nonlinear problems, so that one usually has to be content with using sub-optimal feed- 
back laws. On the other hand, in our study, we consider optimal controllers as determined 
by optimal control methodologies instead of a feedback strategy. 

Optimal control techniques will not provide the real time control that one is ultimately 
interested in, but by systematically computing the best control within specified tolerances, 
with a given objective function, it will be possible to develop strategies (active or passive) 
to control a wide variety of disturbances. For example, to effectively control boundary-layer 
transition due to the interaction of a crossflow vortex and a TS wave using periodic heating 
and cooling, optimal control would allow: (1) a determination of the best objective function 
to use for a given type of control (some are better than others) and (2) provide insight into 
the relationship between the time dependence of the control and the input waves. This 
insight could then be built into a neural network, or other type of self-learning system, to 
allow effective control over a wide range of input parameters. 

1.3 The Current Approach 

The methodology of the current paper is based on defining a control mechanism and 
an objective for control, and then finding, in a systematic and automated manner, controls 
that best meet the objective. In the present setting, an objective or cost functional is 
defined that measures the “distance” between the measured stresses, and their desired 
values along a limited section of the bounding wall and over a specified length of time. 
One may interpret the objective functional as a “sensor,” i.e., the objective functional 
senses how far the flow stresses along the wall are from the corresponding desired values. 
To control the flow, time-dependent injection and suction are imposed along a small orifice 
in the bounding wall. Although the spatial dependence of the suction profile is specified (for 
simplicity), the optimal control methodology determines the time- variation of this profile. 
However, unlike feedback control methodologies wherein the sensed data determines the 
control through a specified feedback law or controller, here the time-dependence of the 
control is the natural result of the minimization of the objective functional. This scenario 
is shown in Fig. 1. We have a sensor that feeds information to a controller that in turn 
feeds information to the actuator. However, in the optimal control setting, the sensor is 
actually an objective functional and the controller is a coupled system of partial differential 
equations that determine the control that does the best job of minimizing the objective 
functional. The present active-control approach is demonstrated for the evolution and 
automated control of spatially growing 2D disturbances in a flat-plate boundary layer. As 
the length of time over which the minimization process is increased, the results obtained by 


6 



wave cancellation are recovered, thus validating the approach. The ultimate goal of this line 
of research is to introduce automated control to external flows over realistic configurations 
such as wings and fuselages, and to devise novel flow control systems. 


2. The Optimization Problem 
2.1. The state equations 

Let 0 denote the flow domain which is the semi-infinite channel or boundary layer 
\x > 0, 0 < y < h\, where h is the location of the upper wall for the channel or the 
truncated freestream distance for the boundary layer. Let , denote its boundary and let 
(0, T) be the time interval of interest. The inflow part of the boundary \x = 0, 0 < y < h\ 
is denoted by , i and the part of the boundary on which control is applied (i.e., along 
which the suction and blowing actuator is placed) by , a which is assumed to be a finite 
connected part of the lower boundary (or wall) \x > 0, y = 0]. Solid walls are denoted 
by , w ; for the channel flow, , w is the lower boundary \x > 0, y = 0] with , a excluded 
and the upper boundary \x > 0, y = h\; for the boundary layer flow, , w is only the lower 
boundary with , a excluded. For the boundary-layer case, the upper boundary \x > 0, 
y = h\, which is not part of , W} is denoted by , e . Controls are only activated over the 
given time interval t £ (To, Ti), where 0 < To < T\ < T. 

The flow field is described by the velocity vector (u, v) and the scalar pressure p and 
is obtained by solving the following momentum and mass conservation equations 


du du du dp d f du\ d f du dv 
dt + dx + dy + dx dx\dx) dy\dy + dx 


0 in (0, T) x ft , (1) 


dv dv dv dp d f du dv 

dt + dx ^ dy + dy dx\dy + dx 


d ( ^dv 
dy V dy 


du dv 
dx dy 


0 


subject to initial and boundary conditions: 


in (0, T) X fi, 


0 in (0, T) xfi, (2) 

(3) 


(u,v)\ t=0 = (uq,vq) in 0 

(4) 

. f ( 91 , 92 ) in (To, Ti) 

(5) 

(u,c) = < 

lr “ [(0,0) in (0,T o ) and (T l5 T) 

u, v) r = (it;, Vi) and (u,u)| r =(0,0) in (0,T) 

(6) 

t du dv 

(7) 

it, v, p) — X base flow, — — , — xO as x — X 00 . 

dx dx 
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Here, the initial velocity vector (uo(x,y),vo(x,y)) and the inflow velocity vector 
(ui(t,y),Vi(t,y)) are assumed given and the base flow is assumed to be Poiseuille flow 
for the channel case and Blasius flow for the boundary-layer case. The above system holds 
for both the channel and Blasius flow cases; in the latter case, the upper boundary is not 
part of , w and the additional boundary condition 


u\r e = Uoo and p-2v— = Poo in (0,T) (8) 

°y v 

J- e 

is imposed, where U a Q and denote the free-stream flow speed and pressure, respectively. 

The control functions g\ (f, x) and g 2 (t, x) which give the rate at which fluid is injected 
or sucked tangentially and perpendicularly, respectively, through , a are to be determined 
as part of the optimization process. In order to make sure that the control remains bounded 
at To, it is required that 


9i\t=T 0 = gio(x) and # 2 |t=T 0 = 92 o(x ) on , a , 


(9) 


where gio(x) and g 2 o{x) are specified functions defined on , „. Commonly, one chooses 
g w {x) = g 20 (x) = 0. 

2.2. The objective functional and the optimization problem 

Assume that , s is a finite, connected part of the lower boundary \x > 0, y = 0] which 
is disjoint from , a and that (T a ,T&) is a time interval such that 0 < T a < T& < T. Then, 
consider the functional 


r>Tb 


»T 5 


J(u,v,p,g!,g 2 ) = — 


f I \T1 ~ Taf d } dt + f I I 

z JT a J r, z JT a J r. 


|t 2 — n | d, dt 


+ 


Pi 


r -Ti 


>T„ 


dgi 

dt 


T \ 9 i | 2 ^ d, dt + 




< -Ti 


>T„ 


dg2 

dt 


+ | <72 1 2 ^ d, dt , 


( 10 ) 


where g\ and g 2 denote the controls and r a (t,x ) and Tf,(t,x) are given functions defined 
on (T a ,Tft) X , s . Note that since , s is part of the lower boundary of the channel or 
boundary-layer wall, t\ = vdu/dy and r 2 = — p + 2 vdv/dy are the shear and normal 
stresses, respectively, exerted by the fluid on the bounding wall along , s and thus r a and 
Tb may be interpreted as given shear and normal stresses, repectively. Then, the boundary 
segment , s can be thought of as a sensor which measures the stresses on the wall. Thus, 
in (10), , s is the part of the boundary , along which one wishes to match the shear and 
normal stresses to the given functions r a and 77 , respectively, and (T a ,T&) is the time 



interval over which this matching is to take place. (There are no diffilculties, other than 
notational, introduced if one wishes to match each component of the stress vector over a 
different boundary segment and/or over a different time interval.) 

The third and fourth terms in (10) are used to limit the size of the control. Indeed, no 
bounds are a priori placed on g\ or g 2 ; their magnitudes are limited by adding a penalty to 
the stress matching functional defined by the first two terms in (10). The particular form 
that these penalty terms take, i.e, the third and fourth terms in (10), is motivated by the 
necessity to limit not only the size of the controls g\ and gp, but also to limit oscillations. 
The constants op, a 2 , (3 1 , and /?2 can be used to adjust the relative importance of the 
terms appearing in the functional (10). 

The (constrained) optimization problem is given as follows: 

Find u,v,p,g\, and g 2 such that the functional J{u, v, p, g \ , g 2 ) given 
in (10) is minimized subject to the requirement that (l)-(7) and (10) are 
satisfied and , for the boundary -layer flow case, (8) is also satified. 


3. The Optimality System 

We first consider, in Sections 3. 1-3.5, the case of a channel flow; the optimality system 
for the boundary layer flow is considered in Section 3.6. 

3.1. The Lagrangian Functional 

The method of Lagrange multipliers is formally used to enforce the constraints (l)-(3) 
and (5). To this end, the Lagrangian functional 
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C(u,v,p,g 1 ,g 2 ,u,v,p,s 1 ,s 2 ) 

r*T b /* 


CX\ 

2 J T a ) T . 

Ti r 


2 


'To Jr a 


Ti - T a \ 2 d, dt + -y- / / |r 2 - T 6 | 2 d, dt 

j T a J r s 

(%■ + l9il 2 ) <*■ * + y /) St +^I 2 ) A * 


^ ^ du du 

T| — + U — + U 


/o 


dp d / du\ d ( du 

o ,/o \dt ' dx dy ^ dx dx\dx) dyydy^dx)) 
T f „ / dv dv dv dp d f du dv\ <9 / du \ \ 

Jq \ dt dx dy dy dx\dy dx J dy \ dy J J 
[ T f du dv\ 

/ / P \d^ + d^) dfldt 

i o Jn \ ox °y j 

a /•To /> i-T /• 

s\(u — gi)d, dt — / / siud,dt— / s\ud,dt 

Jo J r a Jtt J r a 

a /•To i- t-T i- 

s 2 (v — g2)d, dt — / / S 2 V d, dt — / / S 2 V d, dt 

Jo J r a J Ti J r a 


:n) 


is introduced. In (11), it and v are Lagrange multipliers that are used to enforce the x 
and y-components of the momentum equation (1) and (2), respectively, p is a Lagrange 
multiplier that is used to enforce the continuity equation (3), and si and s 2 are Lagrange 
multipliers that are used to enforce the x and y-components of the boundary condition 
(5), respectively. Note that Lagrange multipliers have not been introduced to enforce the 
constraints (4), (6), (7), and (9), so that these conditions must be required of all candidate 
functions it, u, p, yq, and y 2 . 

Through the introduction of Lagrange multipliers, the constrained optimization prob- 
lem is converted into the unconstrained problem: 

Find u, v, p, g\, y 2; u, v, p, si, and s 2 satisfying (4), (6), (7), and (9) 
such that the Lagrangian functional £(u, u,p,yi,y 2 , u, u,p, si, s 2 ) 
given by (11) is rendered stationary. 

In this problem, each argument of the Lagrangian functional is considered to be an 
independent variable so that each may be varied independently. 

The first-order necessary condition that stationary points must satisfy is that the first 
variation of the Lagrangian with respect to each of its arguments vanishes at those points. 
One easily sees that the vanishing of the first variations with respect to the Lagrange 
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multipliers recovers the constraint equations (l)-(3) and (5). Specifically, 


S_£ SC 

8u ’ Sv 

s_c 

Sp 


0 

0 


x- and (/-momentum equations (1) and (2) 


continuity equation (3) 


x and y components of the boundary condition (5) , 

where 5C/5u denotes the first variation of C with respect to u, etc. 

3.2. The Adjoint Equations 

Next, set the first variations of the Lagrangian with respect to the state variables it, 
v, and p equal to zero. These result in the adjoint or co-state equations. Note that, since 
for the channel flow, candidate solutions must satisfy (4), (6), (7), and (9), one has that 


SC_ 8C_ 

’ 8s2 


<5’nU=o = <5"^U=o = 0 on 0, Sg 2 \t=T o =0 on , a , 
Su\r { = <Hr ; = 0; <Hr w = = 0 for(0,T), 

dSu dSv 


Sp , Su, Sv, 


dx ’ dx 

First, consider SC/Sp = 0 which yields 
°2 [ [ Sp (t 2 - r b ) d, + / 


0 as x — > oo for (0, T) 


( 12 ) 


„ dSp ^ dSp . 

u— 1- v — — ) dCldt = 0 

o Jo, \ vy 


JTa J r 3 

for arbitrary variations Sp in the pressure. Applying Gauss’ theorem then yields that 


» t 5 


OL 2 


T 

J- a 


/ ‘T i* / Qfi \ /* T i* 

Sp(r 2 —Tb)d, — / Spi — 1- — ) dCtdt + / / 8p(uni + vn 2 ) d, dt = 0 


to JQ 


\dx dy 


where n\ and n 2 denote the x and y components, respectively, of the outward normal to 
O along , . Choosing variations Sp that vanish on the boundary , but which are arbitrary 
in the interior O of the flow domain yields that 

f)lj f)l) 

- + - = 0 on(0 ,Dxn. (13, 

Now choosing variations Sp that are arbitrary along the boundary , yields that 

0 on(0, T) x , \, s , (0, T a ) x , s ,(T b ,T) x , s 

u ni +vn 2 ={ / dv \ (14) 

- a 2 \-p + 2u— - r b ) on (l a , l b ) X , s , 
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where , \, s denotes the boundary , with , s deleted. Note that in the above derivation of 
(13) and (14), as in the derivations found below, the boundary integrals at infinity do not 
make any contribution due to the last relation in (12). 

Next, consider 8C/8v = 0 where equation (12) has been used to eliminate boundary 
integrals along , p , w and as x — > oo and an integral over 0 at t = 0. First, variations 
8v that vanish at t = 0, t = T, and in a neighborhood of , are chosen, but which are 
otherwise arbitrary. Such a choice implies that all boundary integrals in (12) vanish giving 


dv ^ du ^ dv dv dv dp 

dt + dy ^ dy dx dy dy 


d f du dv\ 

~ U d^\d^ + fa) 


d / < 97 \ 

dy V dy) 


= 0 in (0, T) x O , 


( 15 ) 


where equation (3) is used to effect a simplification. Next, variations that vanish in a 
neighborhood of , , but which are otherwise arbitrary, are chosen to obtain 


v\t=T = 0 in H. 


( 16 ) 


Now, along , , 8v and d8v/dn may be independently selected, provided that (12) is satisfied, 
where d/dn denotes the derivative in the direction of the outward normal to O along , . If 
8v = 0 and d8v/dn varies arbitrarily along , , then 


v 


0 on(0, T) x , \, s , (0, T a ) x , S ,(T 6 ,T) x , 

«2 (~P + ~ T b^j ° n Tb) X , s . 


(17) 


To see this, note that along the inflow, , p n 2 = 0 and d/dn = —d/dx while along the top 
and bottom boundaries n\ = 0, d/dn = ±d/dy } respectively, and, since 8v = 0, d8v / dx = 
0. Note that (14) and (17) agree on the boundary segments where they simultaneously 
apply. Finally, 8v is arbitrarily chosen along , a to obtain 


/ Qy \ 

s 2 = -pn 2 - viunx + vn 2 ) - vl— + — jnx - 2v—n 2 on (0,T) 


(18) 


Next, consider SC/Su = 0. Applying to the resulting equation the same process that 
led to (15)-(18) yields 


du „du „ dv du du dp 

dt dx dx dx dy dx 


d f du 


-v- 


dx \ dx 


d f dii dv 

dy \dy dx 


= 0 in (0,T) x $1. 


( 19 ) 


12 



u\t=T = 0 in O, 


( 20 ) 


0 on(0, T) x , \, s , (0, T a ) x , s , (T&, T) x , 

u = < ( du \ 

«i I - Ta I on (T a ,T b ) x , s , 


( 21 ) 


and 


\ 

5l = -pm - u(uni + un 2 ) - 2//— ni - //( — + — Jn 2 on (0,T)x, a . (22) 


In deriving (21) we have used the assumption that , s is part of the lower boundary of the 
channel so that along , s we have that n 2 = —1. Again, there is no conflict between (14) 
and (21) along boundary segments on which both apply. 

3.3. The Optimality Conditions 

The only first-order necessary conditions left to consider are SC/8g\ = 0 and 8C/8g2 = 
0. (These conditions are usually called the optimality conditions.) Now, since all candidate 
functions g\ and p 2 must statisfy (9), it follows that 8g\ = 0 and <h/ 2 = 0 at t = To. Then, 
take 8C/8g2 = 0 and apply Gauss’ theorem to remove all derivatives from the variation 
8g2- Choosing variations 8g2 that vanish at t = Ti, but which are otherwise arbitrary, and 
using (18) yields 

< 9 2 o 2 1 ( dv\ . . . . 

+ = -fo[ p + 2u d^J on ( T °’ T !) X , a, (23) 


where (17) and the assumption that , a is part of the lower boundary so that, along , a , 
m = 0 and n 2 = — 1 have been used. Now, choosing variations that are arbitrary at t = T\ 
yields that dg2/dt = 0 along , a at t = T\ so that, invoking (9), g2(t,x) satisfies 


g2\t=T 0 = 92o(x) and 


dg2 

dt 


= 0 

t=Ti 


on 


(24) 


Note that, given p and 7, (23) and (24) constitute, at each point x £ , a , a two-point 
boundary value problem in time over the interval (To,Ti). 

In a similar manner, setting SC/8g\ = 0 yields that 


d 2 gi 
dt 2 


+ 9i 


1 8 du 

Ji\di 


on (T 0 , Ti) x , a 


(25) 


9i|i=T 0 = 9io(X 


and 


dgi 

dt 


= 0 

t=Ti 


on 


(26) 
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3.4. Finite Computational Domains 

In the computations, the semi-infinte domain Q (we are still only considering the 
channel flow case) is replaced by a finite domain Qc defined by the introduction of the 
outflow boundary , 0 given by \x = L, 0 < y < h\. Thus, we have that Qc is the rectangle 
[0 < x < L, 0 < y < h\. The outflow does not require the imposition of boundary 
conditions along the outflow boundary , 0 because a buffer zone 27 is attached to the end 
of the physical computational domain, where the governing equations are parabolized in 
this buffer region. 

A similar treatment of the adjoint variables should have required consideration of 
an infinite domain [— oo<x<oo,0<y<h]. If this had been done, the boundary 
conditions (17) and (21) would not have been obtained along the inflow , i. In fact, the 
inflow boundary , i for the state equation is the outflow boundary for the adjoint equations 
and, conversely, the outflow boundary , 0 for the state equation is the inflow boundary for 
the adjoint equations. This is easily seen by comparing the leading inertial terms of the 
state equations with t increasing and the adjoint equations with t decreasing. Now, on 
both , i and , 0 we have that u > 0 and v ss 0 which is why , i is an inflow boundary and 
, o is an outflow boundary for the state. On the other hand, the fact that t is decreasing 
in the adjoint equations implies that now , i is an outflow boundary and , 0 is an inflow 
boundary for those equations. 

Thus, to be consistent with the treatment of the state equations, the adjoint outflow 
, i should be treated in a manner similar to the above treatment of the state outflow , 0 . 
This treatment of the adjoint outflow does not require the imposition of any boundary 
conditions for the adjoint variables along , i. Finally, since , 0 is an inflow boundary for 
the adjoint equations, one has that 

u = 0 and 7 = 0 on (0, T) x , 0 . (27) 

3.5. The Optimality System for Channel Flow 

We now have in hand the full optimality system for channel flow whose solutions 
determine the optimal states, controls, and adjoint states. These are: 

State equations - (l)-(7) 

Co-state equations - (13)-(17), (19)-(21), (27) 

Optimality equations - (23)- (26) 

Since (18) and (22) merely serve to determine the uninteresting Lagrange multipliers 
S 2 and si, they can be ignored. 
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The state equations are driven by the given initial velocity (rto,uo), the given inflow 
velocity (up Vi ), and the controls (<71 , <72 ) - Indeed, the purpose of this study is to determine 
g 1 and g 2 that optimally counteracts instabilities created upstream of , a . The adjoint 
equations are homogeneous except for the boundary condition along , s , the part of the 
boundary along which we are trying to match the stresses. The data in that boundary 
condition is exactly the discrepancy between the desired stresses r a and 77, and the stresses 
7~i = udu/ dy and T2 = — p + 2 udu/ dy along , s , weighted by the factors aq and aq- The 
equations for the controls are driven by the negative of the adjoint stresses along , a , the 
part of the boundary along which we apply the control, weighted by the factors 1 / and 
1 / /?2 - Of course this division into equations for the state, the adjoint state, and the control 
is really obscured by the fact that equations are all intimitately coupled. 

3.6. The Optimality System for Boundary-Layer Flow 

Following a similar process to that used in Sections 3. 1-3.5 for the channel flow case, 
one may derive an optimality system for the boundary layer flow case. The only difference 
is that in the latter case , w denotes only the lower boundary with , a excluded and that 
the additional boundary condition (8) along the upper boundary , e must be accounted 
for. 

With the new interpretation for , W} one can still define the Lagrangian functional 
(11) and use the constraints (12) on allowable variations; however, due to (8), allowable 
variations are further constrained by 


Hr e 


8p — 2 v 


d8v 

dy 


r e 


0 for (0, T) . 


which implies that, along , e , one may not choose the variations in 8p and d8v/dy inde- 
pendently. Considering, simultaneously, variations in p, v, and dv / dy along , e , one can 
show that 

u = 0 on (0,T) x , e (28) 

Then, letting 8v be arbitrary along , e yields 

p + 2z/— + vv = 0 on (0, T) x , e . (29) 

The resulting system for boundary-layer flow now includes (28) and (29) in addition 
to the channel flow system. 
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4. Numerical Experiments 

Here, the optimal control methodology developed in Section 3 is applied to a boundary- 
layer flow having a single instability wave that can be characterized by a discrete frequency 
within the spectrum. As described by Joslin et ah, 21,28 these discrete small-amplitude 
instabilities can be suppressed through wave cancellation (WC) using known information 
about the wave. Hence, the optimal control is “known” for validation of the present 
DNS/optimal control theory numerical approach in which the instability is to be suppressed 
without any a priori knowledge of said instability. 

We note that although we are testing our methodology and code for the special prob- 
lem of two-dimensional TS wave suppression, that these were developed for the fully 
nonlinear Navier-Stokes system and thus are applicable to the case of nonlinear, three- 
dimensional waves. We would expect that more iterations, and perhaps a refined inter- 
ation procedure, would be needed for convergence in a multi-mode case. Also, we have 
successfully suppressed a single large-amplitude TS wave; space limitations have prevented 
us from presenting these results here. 

The formidable coupled system is solved in an iterative manner. First, the Navier- 
Stokes equations are solved for the state variables, i.e., the velocity field (it, v) and pressure 
p with control information (i.e., no control gi = g 2 = 0 for first iteration). Then co-state 
equations are solved for the adjoint or co-state variables (u,v) and p. Then, using these 
adjoint variables, the controls g\ and <72 are then found by solving the optimality equations. 
The procedure is repeated until satisfactory convergence is achieved. 

The nonlinear, unsteady Navier-Stokes equations and linear adjoint Navier-Stokes 
equations are solved by direct numerical simulation (DNS) of disturbances that evolve 
spatially within the boundary layer. The spatial DNS 29,30 approach involves spectral and 
high-order finite-difference methods and a three-stage Runge-Kutta method 31 for time 
advancement. The influence-matrix technique is employed to solve the resulting pressure 
equation. 20,32 Disturbances are forced into the boundary layer by unsteady suction and 
blowing through a slot in the wall. The buffer-domain technique 27 is used for the outflow 
boundary treatment. 

In the present study only normal injection or suction control is allowed, so that we 
set gi = 0 in (5), (3 1 = 0 in the functional (10), and ignore (25) and (26). Also, we only 
match the normal stress along , s so that we choose aq = 0 in the functional (10) and in 
( 21 ). 
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4.1. Computational Parameters 

For the computations, the grid has 401 streamwise and 41 wall-normal points. The 
free stream boundary is located 75$* from the wall, and the streamwise length is 2243* 
which is equal to approximately 8 TS wavelengths. The nondimensional frequency for 
the forced disturbance is F = to / R X 10 6 = 86; the forcing amplitude is Vf = 0.1%. 
The Reynolds number based on the inflow displacement thickness ($*) is R = 900. (The 
boundary segment along which disturbance forcing and control is effected as well as where 
stress matching occurs are located within the unstable region of the linear stability neutral 
curve.) A time-step size corresponding to 320 steps per period T p is chosen for a three- 
stage Runge-Kutta method. Based on the disturbance frequency, a characteristic period 
can be defined as T p = 2tt / u = 81.1781; the resulting time-step size is then At = 0.2537. 

To complete one period of the active-control simulation process, 0.75 minutes on the 
Cray C-90 are required using a single processor. Note, two periods of cost (T a — >■ T& and 
T}, — y T a ) are required to complete one iteration of the DNS/adjoint system. Although in 
general any time interval may be specified for T a — >■ Tj, this study uses integer increments 
of the period (T p ) for simplicity. Hence, T a — > T& = 2 T p would cost 4 T p in computations, 
or roughly 3 min of C-90 time per iteration. Because only a single small-amplitude wave 
(linear) is forced, the above grid is more than adequate; however, a grid refinement was 
performed and produced results equivalent to the results reported here. 

For this study, the disturbance forcing slot , j, the control or actuator orifice , a , 
and the matching or sensor segment , s have equal length 4.48$*. The forcing is centered 
downstream at 389.62$* (the Reynolds number based on the displacement thickness at 
that location is R = 1018.99), the actuator is centered at 403.62$* (R = 1037.13), and the 
sensor is centered at 417.62$* ( R = 1054.97). These separation distances were arbitrarily 
chosen for this demonstration. In practice, the control and matching segments should have 
a minimal separation distance so that the pair can be packaged as a single unit, or bundle, 
for distributed application of many bundles. 

4.2. Results 

All simulations allow the flow field to develop for one period, i.e., from t = 0 — >■ T a = T p 
before control is initiated. In the first series of simulations, the interval during which control 
is applied is arbitrarily chosen to be T a — > T& = 2 T p . Based on aq = = 0, «2 = 1, 

and /?2 = 10, the convergence history for the wall-normal velocity and measured normal 
shear T 2 are shown in Fig. 2. The velocities are obtained at a fixed distance from the 
wall corresponding to 1.18$* and at the fixed time Tj. Convergence is obtained with 4 
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iterations. The results demonstrate that a measure of wave cancellation can be obtained 
from the DNS/control theory system. The wall-normal amplitude of the modified wave 
at R = 1092.5 is 40 percent of the uncontrolled wave; the control without optimizing the 
choice of op, a 2 , (3 1 , and /U has led to a 60 percent decrease in the amplitude of the 
travelling wave. Clearly, Fig. 2 shows that a net reduction of the disturbance energy is 
obtained by energy input due to the control. This results in a delay of transition by-way-of 
a suppression of the instability evolution. 

In the simulation, the control has been applied from T a — >■ T& only; therefore, for 
t > T&, equation (5) indicates that the actuation is discontinued. Figure 3 compares 
velocity profiles at R = 1073.2 for the converged results (Cl) of Fig. 2 with results for one 
period after control t = Tf, -\-T p . The measured disturbance tends toward the uncontrolled 
solution when the actuation is discontinued (as expected); because the control was applied 
for T a — T p — y T}, = 2T p} 2 periods are required after T& before the computed solution in 
the window exactly matches the uncontrolled solution. Clearly, the control only removes 
energy from the system (decreases the wave instability amplitude); the resulting profiles 
retain the expected instability profile shape. 

The effect of varying the window size (T a , T&) is demonstrated. The previous converged 
results (Cl) are shown with converged results (C2) for the extended window (T a = T p — > 
T}, = 3 T p ). The results are identical for the first two periods of time and indicate that 
extending the amount of time for control serves to extend control only. This result also 
indicates that one can solve for the optimal control over a given time interval (T a ,T&) 
by breaking up that inerval and solving for the optimal control over a series of smaller 
subintervals. This approach usually leads to substantial savings in CPU and memory 
costs. An additional insight about the present DNS/control theory is gleamed by the 
increased temporal frame. The resulting optimal control g 2 approaches the desired wave- 
cancellation time-periodic solution as the temporal length (T a ,T&) is increased. This is 
convincing evidence that the present self-contained methodology is valid. 

The instability wave resulting from wave-cancellation (WC) is shown with the control 
(C2) in Figs. 4 and 5. For the present comparison, the amplitude of the actuation for WC 
was adjusted until nearly exact wave cancellation was achieved. Although the DNS/control 
theory did not achieve the same level of energy removal, the similar effect of WC was 
achieved without any a priori knowledge of the instability. Also, note that Fig. 5 shows 
that the optimal control of the control theory has nearly the exact phase characteristics 
as WC and only lacks the necessary amplitude for additional wave cancellation. These 
encouraging results suggest that by the appropriate selection of aq, a 2 , /A, and /U, the 
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optimal control can be made nearly as effective method of instability suppression as exact 
wave cancellation. 

From the wave-cancellation study of Joslin et al., 21 the relationship between ampli- 
tude of the actuator (v a ) with resulting instabilty can be shown in Fig. 6. A similar result 
was shown in the channel flow wave-cancellation study in Biringen. 14 The trend indicates, 
that beginning with a small actuation amplitude, as the actuation level is increased, the 
amount of wave cancellation (energy extraction from the disturbance) increases. At some 
optimal actuation, nearly exact wave cancellation is achieved for the instability wave. As 
the actuation amplitude further increases the resulting instability amplitude increases; this 
was clearly explained in Joslin et al. 21 to occur because in the wave superposition process, 
the actuator wave becomes dominant over the forced wave. At this point, the resulting 
instability undergoes a phase shift corresponding to the phase of the wave generated by the 
actuator. The relationship depicted in Fig. 6 is encouraging for the DNS/optimal control 
theory approach and suggests that a gradient descent type algorithm might further enhance 
the wave suppression capability of the present approach. Namely, an approach for the op- 
timal selection of op, a 2, /A, and /?2 might lead to a more useful theoretical/computational 
tool for flow control. 

To simply demonstrate this concept, Lagrange interpolation (or perhaps extrapola- 
tion) is introduced for /% and /?2 based on imposed values for a 1 and a^'- 


in + 1 _ ^r,2( r l*,2 r g2 1 ) fll, 2 1 ( T 1,2 T l,2) 


on+l- _ 
VI, 2 — 


<T n - I-™ -1 ' 

hi, 2 '1,2 y 


(30) 


where t* 2 are some desired values of the stress components and t ™ 2 are the stress com- 
ponents based on the choice (3™ 2 . Although t* and r 2 * may be equivalent to the target 
values T a and r& in the functional (10), this may lead to significant over/under shoots for 
the interation process. Instead, r* and r 2 * is the incremental decrease, or target value, for 
interpolation to more desirable and /?2 values. To illustrate this process, the (3 2 = 10 
(C2) and /?2 = 11 control results are obtained with the iteration procedure. The measures 
of normal stress are somewhat arbitrarly obtained at some time as measured by the sensor 
or matching segment , s ; the values of the normal stress are given in the Table 1. These 
values are used for a desired normal stress r 2 *, which in this case is 65% of the /?2 = 11 
results. 


Table 1. Normal stress for two values of /?2 - 


th 

normal stress 

10 

11 

9.369 x 10“ 6 
8.814 x 10“ 6 
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Using the results for /U = 10 and (3 2 = 11 in (30) yields the value /U = 16.5 which 
is used in a simulation to obtain a greater degree of instability suppression. The WC 
results and the enhanced optimal control (C3) solution are shown in Figs. 7 and 8. This 
interpolation approach based on relationship of Fig. 6 indicates that optimizing /U has led 
to results very close to WC. The solutions differ somewhat near t = T a and t = T& because 
of the conditions (24) and (26) that serve to control the levels of g\ and < 72 - For all practical 
purposes, the solutions obtained with the present DNS/control theory methodology yield 
the desired flow control features without prior knowledge of the forced instability. 

The adjoint system requires that the velocity field (u } v) obtained from the Navier- 
Stokes equations be known for all time. For the iteration sequence and a modestly course 
grid, 82 Mbytes of disk (or runtime) space are required to store the velocities at all time 
steps and for all grid points. For T a — > T& = 3 T p , 246 Mbytes are necessary for the compu- 
tation. Clearly for three-dimensional problems the control scheme becomes prohibitively 
expensive. Therefore, a secondary goal of this study is to determine if this limitation can 
be elliminated. 

Because the characteristics of the actuator (g 1 and g^) and resulting solutions are 
comparable to WC, some focus should be placed on elliminating the enormous memory 
requirements discussed above. This limitation can easily be removed if the flow-control 
problem involves small- amplitude unsteadiness (or instabilities). The time-dependent co- 
efficients of the adjoint system (15) and (19) reduce to the steady-state solution and no 
additional memory is required over the Navier-Stokes system in terms of coefficients. This 
has been verified by a comparison of a simulation with steady coefficients compared with 
the C2 control case. The results for both cases are identical (as expected). Additionally, 
if the instabilities have small amplitudes, then a linear Navier-Stokes solver can be used 
instead of the full nonlinear solver, which was used in the present study. This linear system 
would be very useful for the design of flow-control systems. However, if the instabilities in 
the flow have sufficient amplitude to interact nonlinearly, then some measure of unsteady 
coefficient behavier is likely required. Depending on the amplitudes, the coefficients saved 
at every time-step may be replaced with storing coefficients every 10 or more time-steps 
thereby reducing the memory requirements by an order of magnitude. This hypothesis 
will require validation in a future study. 

Conclusions 

The coupled Navier-Stokes equations, adjoint Navier-Stokes, and optimality condi- 
tion equations were solved and validated for the flow-control problem of instability wave 
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suppression in a flat plate boundary layer. By solving the above system, optimal con- 
trols were determined that met the objective of minimizing the perturbation normal stress 
along a portion of the bounding wall. As a result, the optimal control was found to be 
an effective means for suppressing two-dimensional, unstable Tollmien-Schlichting travel- 
ling waves. The results indicate that the DNS/control theory solution is comparable to 
the wave-cancellation result but, unlike the latter, requires no a priori knowledge of the 
instability characteristics. 
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C - Optimal Control Theory 


Fig. 1. Schematic of active flow control using optimal control theory. 
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Fig. 3. Disturbance velocity profiles for no control (T=4), control (Cl), and after control 
is used and turned-off (C=4) in flat-plate boundary-layer flow. 











Fig. 6. Disturbance velocity resulting from variations in actuator amplitude from simula- 
tions in Joslin et al. 21 









